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We consider the influence of quenched spatial disorder on phase transitions in classical and quantum sys- 
tems. We show that rare strong disorder fluctuations can have dramatic effects on critical points. In classical 
systems with sufficiently correlated disorder or in quantum systems with overdamped dynamics they can 
completely destroy the sharp phase transition by smearing. This is caused by effects similar to but stronger 
than Griffiths phenomena: True static order can develop on a rare region while the bulk system is still in the 
disordered phase. We discuss the thermodynamic behavior in the vicinity of such a smeared transition using 
optimal fluctuation theory, and we present numerical results for a two-dimensional model system. 
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1 Introduction 

The effects of impurities, dislocations, grain boundaries or other types of quenched disorder on phase 
transitions and critical points have fascinated physicists for more than 3 decades. In the early days, it 
was thought that quenched disorder destroys any critical point, because in the presence of defects, the 
system divides itself up into spatial regions which independently undergo the phase transition at different 
temperatures. At such a smeared phase transition, a unique critical temperature for the entire system 
does not exist, and the singularities in the thermodynamic quantities, which are the characteristic feature 
of a critical point, are rounded (see Ref. [1] and references therein). However, subsequently it became 
clear that generically a phase transition remains sharp in the presence of quenched disorder, at least for 
classical systems with short-range disorder correlations. The fate of a particular clean critical point under 
the influence of impurities is controlled by the Harris criterion [2] : If the correlation length critical exponent 
v fulfills the inequality v > 2/d where d is the spatial dimensionality, the disorder does not affect the 
critical behavior. If the Harris criterion is violated, the generic result in classical systems is a new critical 
point with conventional power law scaling but new exponents which fulfill the Harris criterion. 

In recent years, the problem of quenched disorder and critical points has reattracted a lot of attention be- 
cause new results have challenged the simple classification outlined above. At quantum phase transitions, 
i.e., transitions occurring at zero temperature as a function of a nonthermal control parameter, quenched 
disorder can lead to exotic critical points where the conventional power-law scaling is replaced by ex- 
ponential (activated) scaling (see Ref. [3] for a prototypical example, the random transverse field Ising 
model). Similar effects can also be found in classical systems with linear defects [4]. Recently, even more 
dramatic effects have been proposed to occur at some quantum phase transitions in itinerant electronic 
systems [5]. At these transitions a true static order parameter can develop on spatial regions that are more 
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strongly coupled than the bulk system because of disorder fluctuations. Thus, these rare regions undergo a 
phase transition independently from the rest of the system, leading to a smeared global transition exactly 
in the sense discussed at the beginning of this section. 

In this paper we explore this smeared phase transition scenario in more detail. We show that disorder- 
induced smearing of a critical point is a ubiquitous phenomenon occurring in a wide variety of systems 
ranging from the above-mentioned itinerant electronic systems to classical Ising magnets with planar de- 
fects and to non-equilibrium phase transitions. To illustrate the smearing mechanism, we present computer 
simulations of a two-dimensional (2D) model system. The paper is organized as follows: In Section|2| we 
discuss scaling scenarios for a critical point in the presence of quenched disorder. Section [3]is devoted to 
rare regions and Griffiths singularities. In Section|4]we discuss how these rare regions can destroy a critical 
point, and we present numerical results in Section|5] 

2 Critical point scenarios 

In this section, we discuss possible scaling scenarios for (sharp) critical points in the presence of weak 
quenched disorder. Weak disorder refers to types of defects which introduce spatial variations of the 
coupling strength but no frustration or random external fields. For these systems, three broad classes can 
be distinguished according to the behavior of the disorder under coarse graining. 

The first class consists of clean critical points that fulfill the Harris criterion [2] v > 2/d, At these 
transitions, the disorder decreases under coarse graining, and the system becomes asymptotically homoge- 
neous at large length scales. Consequently, the critical behavior of the dirty system is identical to that of the 
clean system. Technically, this means the disorder is renormalization group irrelevant, and the clean renor- 
malization group fixed point is stable. In this first class the macroscopic observables are self-averaging at 
the critical point, i.e., the relative width of their probability distributions vanishes in the thermodynamic 
limit [6,7]. A prototypical example in this class is the 3D classical Heisenberg model whose clean corre- 
lation length exponent is v w 0.698 (see, e.g., [8]), fulfilling the Harris criterion. 

The other two classes occur if the clean critical point violates the Harris criterion. In the second class, 
the system remains inhomogeneous at all length scales with the relative strength of the inhomogeneities 
approaching a finite value for large length scales. The resulting critical point still displays conventional 
power-law scaling but with new critical exponents which differ from those of the clean system (and fulfill 
the Harris criterion). These transitions are controlled by renormalization group fixed points with finite dis- 
order. Macroscopic observables are not self-averaging, the relative width of their probability distributions 
approaches a size-independent constant [6,7]. An example in this second class is the classical 3D Ising 
model. Its clean correlation length exponent, v « 0.627 (see, e.g. [9]) does not fulfill the Harris criterion. 
Introduction of quenched disorder, e.g., via dilution, thus leads to a new critical point with an exponent of 
v « 0.684 [10]. 

At critical points in the third class, the relative magnitude of the inhomogeneities increases without 
limit under coarse graining. The corresponding renormalization group fixed points are characterized by 
infinite disorder strength. At these infinite-randomness critical points, the power-law scaling is replaced 
by activated (exponential) scaling. The probability distributions of macroscopic variables become very 
broad (even on a logarithmic scale) with the width diverging with system size. Consequently, averages 
are often dominated by rare events, e.g., spatial regions with atypical disorder configurations. This type 
of behavior was first found in the McCoy- Wu model, a 2D Ising model with bond disorder perfectly 
correlated in one dimension [4]. However, it was fully understood only when Fisher [3] solved the ID 
random transverse field Ising model by a version of the Ma-Dasgupta-Hu real space renormalization group 
[11]. Since then, several infinite-randomness critical points have been identified, mainly at quantum phase 
transitions since the disorder, being perfectly correlated in (imaginary) time, has a stronger effect for 
quantum phase transitions than for thermal ones. Examples include ID random quantum spin chains as 
well as ID and 2D random quantum Ising models [12-16]. 
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3 Griffiths singularities 

In the last section, we have discussed the global (average) behavior of the disorder at a critical point. In 
this section, we focus on the effects of rare strong spatial disorder fluctuations. Such fluctuations can lead 
to very interesting non-perturbative effects not only directly at the phase transition but also in its vicinity. 
These effects, which are known as Griffiths phenomena, can be understood as follows: Generically, the 
critical temperature T c of a disordered system is lower than its clean value, T c °. In the temperature interval 
T c < T < T c °, the bulk system is in the disordered phase. However, in an infinite size sample, there is 
an exponentially small, but nonzero probability for finding an arbitrary large region devoid of impurities. 
Such a rare region (or Griffiths island) can develop local order while the bulk system is still disordered. 
The dynamics of a rare region is very slow because flipping it requires a coherent change of the order 
parameter over a large volume. Griffiths [17] showed that the presence of these locally ordered islands 
produces an essential singularity [17, 18] in the free energy in the whole region T c < T < T c °, which is 
now known as the Griffiths region or the Griffiths phase [19]. In generic classical systems the Griffiths 
singularity is weak, and it does not significantly contribute to the thermodynamic observables. To the 
best of our knowledge, the thermodynamic Griffiths singularities in classical systems have not yet been 
observed in experiments. In contrast, the long-time dynamics is dominated by the rare regions. Inside 
the Griffiths phase, the spin autocorrelation function C(t) decays as lnC(t) ~ — (lni) d /( d_1 ) for Ising 
systems [19-23] and as lnC(i) ~ — t 1 / 2 for Heisenberg systems [22,24]. 

In quantum systems or classical systems with perfectly correlated disorder (like the above-mentioned 
McCoy-Wu model), the Griffiths effects are generically stronger. Effectively, the rare regions are now 
extended objects which are infinite in the correlated space or imaginary time directions. This makes their 
dynamics even slower and so increases their effects [25]. In ID and 2D random transverse field Ising 
models, the singularity of the free energy in the Griffiths phase takes a power-law form with nonuniversal 
continuously varying exponents. Several thermodynamic observables including the average susceptibility 
actually diverge in a finite region of the disordered phase [3, 14-16] rather than only at the critical point. 
Similar phenomena have also been found in quantum Ising spin glasses [26-28]. Because quantum Grif- 
fiths singularities are much stronger than classical ones they should be accessible in experiments. Indeed, 
the peculiar behavior of many heavy fermion systems in the vicinity of a magnetic quantum phase tran- 
sition has been interpreted as quantum Griffiths behavior [29] even though the situation is probably more 
complex, at least for Ising symmetry [30], as will be explained in the next section. 

4 Disorder induced smearing of a phase transition 

4. 1 Smearing Mechanism 

In the preceeding sections, we have assumed that the phase transition remains sharp in the presence of 
quenched disorder. We now show that in some systems the rare regions effects can become so strong 
that the transition is completely destroyed by smearing. Consider a single rare region that is locally in 
the ordered phase. In the generic case of uncorrected disorder it is of finite extension. Therefore it 
cannot undergo a true phase transition independently of the bulk system. Instead, it very slowly fluctuates 
leading to the Griffiths effects discussed in section [3] This argument suggests that the phenomenology 
of the transition can completely change, if the system permits rare regions that can undergo a true phase 
transition independently of the bulk, leading to a static order parameter on the rare regions. It was recently 
found [5, 31-33] that this does indeed happen, and that it leads to a smeared global phase transition. 

The behavior of such a smeared transition is markedly different from that of a conventional (sharp) con- 
tinuous phase transition. At a conventional transition, a non-zero order parameter develops as a collective 
effect of the entire system which is signified by a diverging correlation length of the order parameter fluctu- 
ations at the critical point and accompanying singularities in the thermodynamic observables. In contrast, 
at a smeared transition the system divides itself up into spatial regions which independently undergo the 
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transition at different values of the control parameter. Once true static order has developed on some of the 
rare regions their order parameters can be aligned by an infinitesimally small interaction or external field. 
Therefore, the global order develops very inhomogeneously over a range of control parameter values. The 
correlation length remains finite, and the singularities of the thermodynamic observables are rounded. To 
show the ubiquity of this smeared transition scenario we now discuss three different examples. 



4.2 Magnetic quantum phase transitions of itinerant electrons 

The first class of systems displaying smeared transitions are itinerant quantum magnets with Ising order 
parameter symmetry. For definiteness we consider the antiferromagnetic quantum phase transition. The 
Landau-Ginzburg-Wilson free energy functional of the clean transition reads [34, 35] 

S = J dx dy m(x)T(x,y) m(y) + u J dx m 4 (x) . (1) 

Here m is the staggered magnetization, x = (x, r) comprises position x and imaginary time t, and 
J dx = J dx J^ T dr. Y(x,y) is the bare two-point vertex, whose Fourier transform is T(q, ui n ) = 
(t + q 2 + \u! n \). Here t = (g — g c ) / g c is the distance of the coupling constant g from the (clean) critical 
point. The dynamical part of Y is proportional to |a> n | reflecting the overdamping of the dynamics due 
to the coupling of the order parameter to fermionic particle-hole excitations (undamped dynamics leads 
to w 2 ). Quenched disorder is introduced by making t a random function of position, t — > t + Jt(x). 
We consider a Poisson (dilution) type of disorder, St(x) = everywhere except on randomly distributed 
finite-size islands (impurities) of spatial density p where 6t(x) — W > 0. 

The rare regions in this system are large spatial regions devoid of these impurities. If t < 0, some of 
the rare regions are locally in the ordered phase. The crucial difference of the itinerant system considered 
here and the transverse field Ising models of section|3]is the overdamped dynamics. The linear frequency 
dependence in Y is equivalent to a long-range interaction in imaginary time of the form (r — r')~ 2 . Each 
rare region is thus equivalent to an ID Ising model with a 1/r 2 interaction. This model is known to have 
a phase transition [36]. Thus, true static order can develop on those rare regions which are locally in the 
ordered phase, leading to a smeared global phase transition as described above [5]. Physically, this means 
the rare regions become static because the damping prevents them from tunneling, as was independently 
shown in Ref. [30]. 



4.3 Classical magnets with extended defects 

A second class of systems which can show smeared transitions are classical magnets with extended impu- 
rities. As an example, consider a 3D classical Ising model with planar defects, given by the Hamiltonian 

H = — JiSi t j,kSi-\-ij,k — J(Sij t kSij-\-i f k + Sij^Sij^+x) ■ (2) 

i=l,...,L ± i=l L ± 

j , fc = 1 , . . . , L j , k= 1 , . . . , L q 

Here, Syfc = ±1 are classical Ising spins on the sites of a cubic lattice. In the clean system all interactions 
are identical and have the value J. The defects are modelled via 'weak' bonds randomly distributed in one 
dimension, the uncorrelated dimension. The bonds in the remaining two (correlated) dimensions remain 
equal to J. Lj_ and Lc are the system lengths in the uncorrelated and correlated directions, respectively. Ji 
is the random coupling constant in the uncorrelated direction, drawn from a binary probability distribution: 



P[Ji] = (l-p)6(Ji-J)+p8(Ji-cJ) (3) 

where p and c are constants between and 1 . Thus, the Ising magnet has planar impurities of concentration 
p and relative strength c of the weak bonds. 
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Rare regions in this system consist of large blocks of planes devoid of weak bonds. They are infinite in 
two dimensions (the correlated dimensions) but finite in the uncorrelated direction. Since a 2D Ising model 
has an ordered phase, each rare region can independently undergo the magnetic phase transition. This 
leads to a smeared global phase transition in the 3D Ising model with planar defects [31, 32]. Physically, 
the reason for the smearing is different from the itinerant magnets discussed above. Here, an isolated rare 
region can develop a static order parameter because its dimensionality is above the lower critical dimension 
d~ = 1 for the classical Ising model. This also shows that a continuous order parameter symmetry 
suppresses the tendency towards smearing. In classical magnets with XY or Heisenberg symmetry and 
short-range interactions, the defects must be at least three-dimensional for smearing to occur since the 
lower critical dimension for these systems is d~ =2. 

4.4 Nonequilibrium transitions 

A third class of smeared phase transitions are non-equilibrium phase transitions with extended spatial de- 
fects. Consider, e.g., the contact process [37], a prototypical system in the directed percolation universality 
class [38]. It is defined on a rf-dimensional hypercubic lattice. Each site can be vacant or active, i.e, occu- 
pied by a particle. During the time evolution, particles are created at vacant sites at a rate An/ (2d) where 
n is the number of active nearest neighbor sites and the 'birth rate' A is the control parameter. Particles are 
annihilated at unit rate. For small A, annihilation dominates, and the absorbing state without any particles 
is the only steady state. For large A there is a steady state with finite particle density (active phase). Both 
phases are separated by a nonequilibrium phase transition in the directed percolation universality class at 
A = A°. We introduce quenched spatial disorder by making the birth rate A a random function of the lattice 
site. As in the preceeding section, extended impurities can be described by disorder perfectly correlated in 
dc dimensions, but uncorrelated in the remaining d R = d — dc dimensions. A is thus a function of 
which is the projection of the position vector r on the uncorrelated directions. For definiteness, we assume 
that the A(r^) have a binary probability distribution 



where p and c are constants between and 1 . In other words, there are extended impurities of density p 
where the birth rate A is reduced by a factor c. 

Rare regions in this disordered contact process are infinite in the correlated dimensions dc > 1 but 
finite in the d R random directions. The directed percolation universality class permits a phase transition 
for all dimensions d > 1. Therefore, each rare region can undergo a phase transition independently of the 
bulk system. Following the general arguments above, we thus conclude that the global phase transition of 
a system in the directed percolation universality class is smeared by extended spatial defects [33]. 

4.5 Extremal statistics theory 

In this section we derive the leading behavior in the tail of a smeared phase transition, i.e., in the parameter 
region where a few isolated rare regions have already ordered, but their distance is so large that they 
are effectively independent. For definiteness we consider the first example, a magnetic quantum phase 
transition of itinerant electrons. Adaption to the other discussed systems is straight forward. The theory is 
similar to that of Lifshitz [39] and others for the band tails in a disordered semiconductor. 

The probability w for finding a region of linear size L RR devoid of any impurities is given by w <~ 
exp(—pL RR ) (up to pre-exponential factors). Such a rare region develops static order at a distance 
tc(LRn) < from the clean critical point. Finite size scaling yields \t c (L RR )\ ~ L RR where <fi is the 
finite-size scaling shift exponent of the clean system [40]. Thus, the probability for finding a rare region 
which becomes critical at t c is given by 



P[\{v R )] = (l-p) 5[\{v R ) - A] + p 6 [A(r«) - cA] 



(4) 



w{t c ) ~ exp(-B |i c |- dAi ) fort^O 



(5) 
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The total order parameter m is obtained by integrating over all rare regions which are ordered at t, i.e., all 
rare regions having t c > t. This leads to an exponential tail of m as a function of the distance t from the 
clean critical point: 



log(m) ~ — B\t 



-d/4> 



for t -» 



(6) 



Note that the functional dependence on t of the order parameter on a given island is of power-law type and 
thus only influences the pre-exponential factors. 

At finite temperatures, the static order on the rare regions is destroyed, and a finite interaction of the 
order of the temperature is necessary to align them. This means a sharp phase transition is recovered. To 
estimate the transition temperature we note that the interaction between two rare regions depends exponen- 
tially on their spatial distance r which itself, according to l|5}, depends exponentially on t. This leads to a 
double-exponential dependence of the critical temperature on t 



log(-alogT c ) 



1*1 



-d/<f> 



(7) 



We now turn to finite-size effects at zero temperature. The total order parameter is the sum of contribu- 
tions of many independent islands. Therefore, finite size-effects in a macroscopic sample are governed by 
the central limit theorem. However, for t — > 0—, very large and thus very rare regions are responsible for 
the order parameter. The number N of rare regions which start to order at t c in a sample of size L behaves 
like N ~ L d w(t c ). When N becomes of order one, strong sample-to-sample fluctuations arise. Using eq. 
for w(t c ), we find that sample-to-sample fluctuations become important at 



\t L \ 



(logi) 



-<t>/d 



(8) 



The spatial magnetization distribution in the tail of the smeared transition is very inhomogeneous. On 
the ordered rare regions, the local order parameter m(r) is comparable to that of the clean system. Away 
from these islands, it decays exponentially with distance. The probability distribution P[logm(r)] is there- 
fore very broad, ranging from logm(r) = 0(1) on the largest islands to logm(r) — > — oo on sites very 
far away from an ordered island. The typical local order parameter m ty p can be estimated from the typical 
distance of any point to the nearest ordered island. From (|5} we obtain 



r typ ~ exp(B \t\- d ^/d) . 
At this distance from an ordered island, the local order parameter has decayed to 



exp 



-CeMB \t\- d/ */d) 



(9) 



(10) 



where £o is the bulk correlation length (which is finite and changes slowly throughout the tail region of 
the smeared transition) and C is constant. A comparison with l|6) gives the relation between m typ and the 
thermodynamic (average) order parameter m, 



logm 



typ\ 



-l/d 



(ID 



Thus, mtyp decays exponentially with m indicating an extremely broad local order parameter distribu- 
tion. In order to determine the functional form of the local order parameter distribution, first consider 
a situation with just a single ordered island at the origin of the coordinate system. For large distances 
r the magnetization falls off exponentially like m(r) = mo e~ r ^°. The probability distribution of 
x — log[m(r)] = log mo — r/£o can be calculated from 



P(\x\ 



dN 


dN 


dr 


dx 


dr 


dx 



dN 
dr 



(12) 
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where dN is the number of sites at a distance from the origin between r and r + dr or, equivalently, 
having a logarithm of the local magnetization between x and x + dx. For large distances we have \x\ ~ r. 
Therefore, the probability distribution of log m generated by a single ordered island takes the form 

P[log(m)] - |log(m)| d_1 (form<l). (13) 

In the tail region of the smeared transition, the system consists of a few ordered islands whose distance is 
large compared to £q. The probability distribution of log[m(r)] thus takes the form (1131 with a lower cutoff 
corresponding to the typical island-island distance and an upper cutoff corresponding to a distance £o from 
an ordered island. 



5 Numerical results 

In this section, we illustrate the disorder-induced smearing of a phase transition by showing numerical data 
for a classical Ising model with two spatial and one time-like dimensions. The disorder is totally correlated 
in the time-like direction. The interactions are short-ranged in space but infinite-ranged in the time-like 
direction. This simplification retains the crucial property of static order on the rare regions but permits 
system sizes large enough to study exponentially rare events. The Hamiltonian reads 

H = — / j Sx,T»Sy,T' T / j ^x^x.tSx.t' • (14) 
Lit Lit 

(x,y),T,T x,r,r 

Here x and r are the space and time-like coordinates, respectively. L T is the system size in time direction 
and (x, y) denotes pairs of nearest neighbors. J x is a quenched binary random variable with the distribution 
P(J) = (1 — p) S( J — 1) + p S(J). In the following, we fix the impurity concentration at p = 0.2, unless 
otherwise noted. In this classical model, L T takes the role of the inverse temperature in the corresponding 
quantum system and the classical temperature takes the role of the coupling constant g. Because the 
interaction is infinite-ranged in time, the time-like dimension can be treated in mean-field theory. For 
L T — ► csc, this leads to a set of coupled mean-field equations (one for each x) 

m x = tanh/3 [J x m x + m y + h] , (15) 
y(x) 

where h = 1CP 8 is a small symmetry-breaking magnetic field. Eqs. (I15> are solved numerically in a 
self-consistency cycle. 

The left panel of Fig.^shows the total magnetization and the susceptibility (corresponding to the inverse 
energy gap of the quantum system) as functions of temperature for linear size L = 100. The data are 
averages over 200 disorder realizations [41]. At a first glance these data suggest a sharp phase transition at 
T « 4.88. However, a closer investigation (right panel of Fig.[0 shows that the singularities are rounded. 
If this rounding was a conventional finite-size effect the magnetization curve should become sharper with 
increasing L and the susceptibility peak should diverge. This is not the case here. Instead, the transition is 
intrinsically smeared for L — » oo. 

For comparison with the analytical results, Fig.|2]shows the logarithm of the average magnetization as a 
function of 1/ (T° — T) for four different impurity concentrations. Here, T c ° = 5 is the critical temperature 
of the clean system (p — 0). All data sets follow eq. l|6) over several orders of magnitude in m with the 
expected shift exponent of <p = 2. The slope of each curve gives the corresponding prefactor B in eq. 0. 
In our lattice model, B should be proportional to — ln(l — p). The inset of Fig.|2]shows that this relation 
is fulfilled in good approximation. 

In Fig.[5]we demonstrate the extreme inhomogeneity of the system in the tail of the smeared transition. 
It shows the local magnetization of a single sample of linear size L — 400 for T — 4.8875 as a function of 
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Fig. 1 Magnetization and susceptibility in the vicinity of the smeared transition. Left: Overview of the behavior for 
L = 100 and p = 0.2. There is a seeming transition at T ~ 4.88. Right: The tail region of the smeared transition. The 
agreement of the results of different systems sizes shows that the observed rounding is not a finite-size effect (from 
reference [5]). 
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Fig. 2 Log(m) as a function of the distance from 
the clean critical point for four different impurity 
concentrations. The solid lines are fits of the aver- 
age magnetization to eq. {5J with = 2. The inset 
shows the dependence of the slope of the curves on 
-ln(l-p). 



the spatial coordinates. Clearly, a sizable magnetization only exists on a few isolated islands; in-between 
it is vanishingly small. Figure |4] shows the probability distribution of the local magnetization for several 
temperatures, calculated from a single sample of linear size L = 2000. In the tail of the smeared transition, 
i.e., for temperatures close to T c ° = 5, the distribution becomes extremely broad, even on a logarithmic 
scale. A more detailed investigation [31,42] shows that this distribution indeed follows eqs. il Q and d!2i . 



6 Summary and conclusions 

To summarize, we have shown that in several disordered classical and quantum systems, isolated rare 
spatial regions can undergo a true phase transition independently of the bulk system. As a consequence, 
different parts of the system undergo the phase transition at different control parameter values, i.e, the 
global phase transition is smeared by the disorder. For classical phase transitions, this scenario requires 
extended defects whose dimensionality is above the lower critical dimension for the universality class in 
question. Specifically, for classical Ising and Heisenberg magnets the (equilibrium) phase transition is 
smeared if the dimensionalities of the defects are at least two and three, respectively. Non-equilibrium 
phase transitions in the directed-percolation universality class are smeared if the defects are at least linear. 
Quantum phase transitions display a greater tendency towards smearing because quenched disorder is 
always perfectly correlated in imaginary time direction (which has to be taken into account at quantum 
phase transitions). At magnetic quantum phase transitions in itinerant electron systems, the rare region 
effects are further enhanced by Landau damping. This leads to smearing even for point-like impurities. In 
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Fig. 4 Probability distribution of the local magne- 
tization of a single sample of linear size L = 2000 
in the random directions for several temperatures. 



this final section we discuss questions of universality, the relation to Griffiths phenomena, and the Harris 
criterion. We also discuss favorable conditions for observing the smearing in experiments or simulations. 

The origins of the disorder-induced smearing of the phase transition described in this paper and of 
Griffiths phenomena are very similar, both are caused by rare large spatial regions which are locally in the 
ordered phase. The difference lies in the dynamics of the rare regions. In the case of Griffiths effects, the 
order parameter on a rare region fluctuates slowly leading to the singularities [17] discussed in Sec. [3] In 
contrast, smearing occurs if the rare regions actually develop true static order. 

A second remark deals with the relation of the disorder-induced smearing and the Harris criterion. We 
emphasize that the phase transition can be smeared by planar defects even if the corresponding clean critical 
point fulfills the Harris criterion and appears to be stable. The reason is that the Harris criterion assumes 
a homogeneous transition and studies the behavior of the coarse-grained (root-mean-square) disorder at 
large length scales. However, the formation of static order on an isolated finite-size rare region is a non- 
perturbative finite-length scale effect in the tail of the disorder distribution. This type of effects is not 
covered by the Harris criterion. 

Let us also note that the functional dependence of the magnetization and other observables on the 
temperature is not universal, it depends on details of the disorder distribution. Therefore, only the presence 
or absence of smearing is universal in the sense of critical phenomena (i.e., depending on dimensionality 
and symmetry only) while the thermodynamic relations are non-universal. 
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Critical points and quenched disorder: Critical points and quenched disorder 



We now turn to favorable conditions for the observation of the smearing in experiments or simulations. 
A significant number of large rare regions will only exist if the defect concentration is small. On the other 
hand, the impurities have to be sufficiently strong so that the bulk system is still far away from criticality 
when the first rare regions start to order. Thus, the most favorable type of disorder for the observation of 
the smearing is a small density of strong impurities. 
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